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Abstract 

Magnetic excitations in cycloidal magnet LiCu202 are explored using THz absorption spec- 
troscopy in magnetic fields up to 12 T. Below the spin ordering temperature we observe eight 
optically active transitions in the spin system of LiCu202 in the energy range from 4 to 30 cm^^. 
In magnetic held the number of modes increases and the electric polarization flop is seen as a 
change in magnetic field dependence of mode energies. The polarization dependence of two of the 
modes in zero magnetic field fits the selection rules for the cycloid tilted by = 41 it 1° from the 
be plane. For the remaining six modes electric and magnetic dipole approximations cannot explain 
the observed polarization dependence. We do not see the electromagnon in the explored energy 
range although there is evidence that it could exist below 4cm^^. 

PACS numbers: 75.30.Ds, 76.50.+g, 75.25. +z, 75.50.-y, 78.30.-j, 78.20. Ci 
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The coupling between magnetic and electric orders in multiferroic materials makes them 
attractive for technological applications. An especially interesting class of multiferroics (MF) 
are frustrated magnets, where the charge order is driven by the magnetic order- and even 
moderate magnetic fields are sufficient to select between different magnetic ground states 
and change the charge order due to magnetoelectric (ME) coupling. The ME interaction 
transfers part of the electric dipole moment of the charge-ordered state to the spin waves 
enhancing their ability to absorb electromagnetic waves, typically in the THz range.-i^ A 
priori, the knowledge of the strength of the static ME coupling for a given material does 
not tell us the strength of the dynamic ME coupling. Therefore the exploration of THz 
properties is of interest to fully exploit MF materials. 

The spin cycloid is an example of a frustrated magnetic structure where ME coupling 
is allowed. In the spin cycloid the electric polarization is induced either through the spin 
current^ or the inverse Dzyaloshinskii-Moriya (DM) interaction mechanism.-^^ In an one- 
dimensional chain of spins along the y axis with spins rotating in the yz plane the electric 
polarization P = (0, 0, P) is perpendicular to the cycloidal ordering vector Q = (0, Q, 0) and 
perpendicular to the normal of the spin rotation plane, P oc Q x (Sj x Sj+i), see Fig. [lb. 
The orientation of the spin rotation plane is determined by anisotropic interactions, like 
easy plane anisotropy or DM interaction. Several multiferroic materials have a cycloidal 
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spin order: TbMnOg, DyMnOg Qffl, NiaVsOg 
CuO Among them, LiCuaOa was the first low spin, S = 1/2, system to show electric 
polarization in the cycloidal phase. Although LiCuaOa seems to be a good model system 
with prevailing one-dimensionality of its magnetism, the connection between the electric 
polarization and the underlying magnetic structure is not well understood. 

The structure of LiCuaOa, shown in Fig.[T^, contains Cu-0 chains extending along the 
crystallographic b axis. The chains form a zig-zag Cu^^ S = 1/2 spin ladder with the 
ladder plane tilted out from the be plane.— The ladders in the planes, stacked along the 
c axis, have an alternating tilt angle ±38.4° leading to two inequivalent ladders per unit 
cell. The first magnetic phase transition is at 24.6 K into the sinusoidal spin structure.— 
An incommensurate cycloidal spin order with Q || b, driven by frustrated isotropic exchange 
interactions inside the ladders, appears in the second phase transition^^ii^ii^ below 23 K 
together with spontaneous polarization^i P || c. The knowledge about the orientation of the 
cycloid plane is controversial. Initially an ab plane cycloid was considered in the neutron 
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FIG. 1: (color online) (a) Crystal structure of LiCu202- Cu^"*" spins, joined by thick black lines, 
form a zig-zag spin ladder, (b) Cycloidal order of spins Sj positioned along the y axis, shown with 
thick solid arrows. Small arrows show the change 5Sj to visualize the spin wave dynamics. The uoy 
and tOz modes change Sj in x direction and are drawn by blue arrows (solid line envelope) and red 
arrows (dashed line envelope), respectively. Green arrows in the yz plane show the i^x mode. 

scattering study.— However, the direction of P together with the inverse DM mechanism 
dictates a he plane cycloid.— A later study by Seki et al— confirmed the be plane cycloid and 
P II c. However, they found deviations from the be cycloid in accord with the NMR study.— 
The presence of ME coupling in LiCu202 is demonstrated by the flop of polarization from 
P II c to P II a if magnetic field Bq > 2 T is applied along b axis.^A Such a flop can be explained 
by invoking a ME coupling between two spin chains in neighboring ladder planes.— Also, 
resonant soft x-ray magnetic scattering shows evidence that the inter-chain spin coupling 
along the c direction is essential for inducing polarization.— 

Modes at wavevectors k = and k = ±Q couple to long-wavelength electromagnetic 
radiation within the linear spin wave model of a spin cycloid. They correspond, as shown in 
Fig.[T)D, to the rotations of spins about three axes by small angles^. The rotation about the 
X axis (k = 0) costs no energy because the spins rotate within the easy plane, hence Ux = 0. 
Rotations about two other axes have finite energies, Uy = uj^, the magnitude determined by 
the size of the easy plane anisotropy. The motion of spins induced by magnetic component 
Bi of THz radiation is described by SSi/St oc Bi x S^. Accordingly, all three modes are 
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FIG. 2: (color online) Light absorption in LiCu202 (sample L3) below 180 cm~^ measured in 
El ||a(b) (dashed red line) and Ei ||c (solid black line) polarizations. T dependence of absorption, 
El II c, is shown by differential absorption spectra in (b) and in the inset (c) . 

magnetic dipole active, Bi excites Uy and excites with equal intensity in zero field; 
B'^ excites although this mode has zero energy. In the presence of dynamic ME coupling 
the Uy mode is excited by the electric component of THz radiation^ and is referred to as 
the electromagnon.— Thus, the orientation of the spin cycloid can be determined from the 
polarization dependence of THz absorption. 

In this Communication we study magnetic field and polarization dependence of the THz 
absorption in the cycloidal phase of LiCu202. Our aim is to find the electromagnon, the 
orientation of the spin cycloid, and changes in the spin excitation spectrum upon the po- 
larization flop in 2T fleld. Previously only one THz absorption study on LiCu202 has been 
reported^^, the experiment being limited to Bq || c, where a peak at 12 cm~^ with increasing 
energy in the fleld was found. 

We use three single crystal samples from the same batch as in Ref. ll]. Sample LI has ab 
plane area of 6.16 mm^ and thickness of 0.75 mm in c direction, L2 and L3 have ac plane areas 
1.44 mm^ and 6.88 mm^ with thicknesses 1.75 mm and 4.3 mm in b direction, respectively. 
Since the crystals are twinned in the ab plane, with a twin size of the order of 10 /im, we 
use a{b) notation when referring to direction-sensitive physical properties. THz absorption 
spectra are recorded using a Martin-Puplett spectrometer with a 0.3 K bolometer and a 
rotatable polarizer in front of the sample.— 
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FIG. 3: (color online) Polarization dependence of absorption in LiCu202 below 35 cm~^, T = 3.5 K, 
Bq = OT. (a) Differential absorption spectra are shown by curves and (b) absorption line areas 
from Lorentzian fits are given by symbols for following polarizations: black line and circles for 
El II a(b),Bi II b(a) is an average of spectra measured on samples LI and L2; red line and down 
triangles for Ei || a(b),Bi || c, sample L2; blue dashed line and up triangles for Ei || c, Bi || a(b), 
sample L2. The letters from A to F in (a) label lines, see TablelH 

The absolute spectra in Fig. [2^ show several strong phonon absorption lines. The differ- 
ential spectrum in Fig. [2b shows that only the lines below 40 cm~^ that disappear above 
the cycloidal ordering temperature T = 23K are spin wave excitations. In the rest of the 
paper we use the T spectrum at 30 K as a reference spectrum by subtracting it from the 
spectra measured in the given field and at lower T. 

Two peaks are emerging already in the 24 K spectrum in the sinusoidal phase and are 
easily seen at 23 K (Fig.lJt). As T drops the peak energies go up and eventually at least 
six peaks can be seen in the 4K spectrum. Higher resolution spectra show that the two 
strongest peaks at 12 and 24 cm^^ are doublets, split by 0.8 cm~^, Fig.O This increases the 
total number of peaks to eight (Table[T]). As Bq is applied. Fig. HI lines shift, change their 
intensity, and some split. In fields stronger than 5T new modes, G and H, emerge into 
the spectral window from the low energy side, Fig.H]:. The expected number of optically 
active modes within the spin wave approximation is two, ujy and uJz- There are two spin 
chains in the unit cell and the coupling between the chains in c direction may split C and E 
modes into C^, and E^, E^. In addition, if the cycloid is not circular, ellipticity activates 
otherwise optically silent modes at multiples of the ordering vector, ±nQ.— However, the 
electric dipole active isotropic exchange modulation mechanism that couples light to the 
spin wave at Brillouin zone boundary^ is apparently not active because optically active spin 
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TABLE I: Spin excitations observed in the THz absorption spectra of LiCu202 in the order of 
increasing energies at T = 3.5 K and Bq = OT. Best fit Lorentzian parameters - center energy Hloq 
(cm~^), hnewidth at half maximum 7 (cm^^) and hne area (cm^^) are hsted for three different 
orientations of Ei and Bi. 



A 
B 

D 

El 

E2 

F 



Ei||a(6) 
Bill 6(a) 



hiOQ 7 A 



6.6 0.3 0.5 

11.4 0.3 0.2 
12.1 0.5 1.5 
19.1 0.8 0.3 
23.8 0.3 0.1 

24.5 0.7 1.4 

28.6 1.2 0.2 



El II c 

Bi \\a{b) 



hujQ 7 A 



6.6 0.4 1.0 

8.3 0.2 0.1 

11.4 0.4 1.0 
12.0 0.5 1.7 
19.2 0.7 0.4 
23.9 0.5 0.7 

24.5 0.6 1.5 
28.4 1.2 0.3 



Ei||a(6) 
Bi II c 



fiujQ 7 A 



6.6 0.3 1.1 

8.3 0.2 0.1 

11.4 0.4 1.8 

12.0 0.5 1.3 

19.2 0.8 0.6 

23.9 0.5 2.0 

24.4 0.6 1.1 

28.4 0.7 0.6 



modes in LiCu202 are below 30cm~ while the zone boundary energy of the spin waves in 
the h direction is more than 60 cm"^, Ref. lisl. 

The polarization dependence of THz absorption spectra is presented in Fig. [3] and fit 
results of the absorption lines are given in Tabled Our measurement rules out both ah and 
he plane spin cycloids within the selection rules discussed above and outlined in Fig.[T]D. For 
the LiCu202 crystal structure a model where there are two cycloids corresponding to two 
zig-zag chains, tilted out from the he plane by an angle ±6', seems natural. The absorption 
line area for the Bi || c polarization becomes Ac = 4ylocos^^^ with Aq standing for the line 
area when Bi is in the cycloid plane. The factor 4 is two cycloids per unit cell times two 
twin domains. For Bi in the ah plane the line area is Aab = 2Ao(sin^ ^ + 1) where the first 
term in brackets is from Bi ± b domains and the second from Bi || b domains. Lines 
and satisfy this model and give 6' = 41 ± 1° for the cycloid tilt angle which is close to 
crystallographic tilt angle 38.4° of the zig-zag chains. But the tilted cycloid model does not 
explain the polarization dependence of other lines. Moreover, the magnetic and the electric 
dipole approximations do not hold for lines A, and E^. They are weakest when light 
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FIG. 4: (color online) Bq dependence of THz absorption in LiCu202 at 3.5 K. Differential absorp- 
tion spectra (a) Bq || c and Ei || a(b), Bi || b(a), sample LI ; (b) Bq || b(a), Ei || c, Bi || a(b), sample 
L2. Panel (c) has the results of Lorentzian fits where the area of the symbol is proportional to the 
absorption line area A; black empty circles Bq || c and filled red circles Bq || b(a). Black solid and 
dash-dot lines are eye guides, except line H which is a linear fit of line position as a function of Bq. 

propagates parallel to the c axis and stronger when light propagates perpendicular to the 
c axis. In the dipole approximation the absorption line area depends on the orientation of 
light polarization with respect to the crystal axes and not on the orientation of the k vector 
of light. A possible explanation of this k dependence is the circular dichroism or electric 
quadrupolar transitions, but this needs further theoretical study. 

Two cases of magnetic resonance in helical structures have been treated theoretically in 
the spin wave approximation, Bq perpendicular to the spin rotation plane and Bq in the 
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plane of spin rotation.— i2Ii2Si2S These results cannot be directly adapted to our experiment 
since the magnetic field would be perpendicular to the cycloidal plane only for an ab or 
be plane circular cycloid in LiCu202. Given the complexity and quantum nature of the 
problem, the theoretical description of spin wave modes in LiCu202 is not within the scope 
of the current paper. The number of modes is larger in Fig|3)D than in Fig.H^ because two 
crystal domains, Bq || b and Bq -Lb, contribute to the absorption; in Bq || c twin domains 
are equal with respect to the field. Some changes are evident in the THz absorption spectra 
around 2 T, Fig.HjD. This is due to reordering of spins when P fiops from c to a axis. Line A 
splits between 1 and 2T and lines and acquire dispersing modes C|| and C|, Fig.Hb. 
We assign the C|| and C| modes together with new modes Ay and Dy to domains where 
Bq II b since this is the field direction where P is fiopped. There is some evidence, Fig.Ht, 
that A|| and G modes anti-cross and therefore the G mode may belong to the Bq || b domains. 

Electromagnons contribute to the zero frequency dielectric constant as Ae = e(0) — £00 = 
J2j{i^p/i^o)'^ , where the sum runs over all spin wave modes observed below 35cm~^ in the 
El II a THz absorption spectrum, ujq and cj^ are the j-th mode resonance energy and plasma 
frequency, both in wavenumber units (cj^ is related to line area by nuj^ = V An with refraction 
index n = 3.36 for Ei || a); 600 contains the contribution from all other modes above 35cm~^ 
and is taken T-independent below 30K.^° Assuming that all the oscillator strength is due 
to electric dipole activity we can set an upper limit to Asa = 0.024 using Ei || a(b), Bi || c 
data from TablelB A measurement at 28 kHz gave Asa = 0.1, Ref.llll. which is four times 
larger. However, the Bq dependence shows (Fig.H]:) that there is an additional mode (H) 
which shifts into our measurement window at linear rate 0.522 ± 0.005 cm^^T^^ with a zero- 
field intercept at 1.03 ± 0.05 cm^^. Assuming the same oscillator strength for the zero-field 
resonance as is above 5 T we get that the contribution of this 1 cm~^ mode to the dielectric 
constant is Asa = 0.46. The actual number might be smaller because it is not known if, first. 
Up = const, and second, whether c^o oc Bq persists down to zero field. A mode at lcm~^, 
when extrapolated to T, has been detected with an electron spin resonance technique 
although with nonlinear Bq dependence.— 

In conclusion, the analysis of the THz absorption spectra within the linear spin wave 
theory shows that the spin order is not ab nor be plane cycloid. The cycloid tilt angle 41 ±1° 
out of the be plane is inferred from the polarization dependence of two modes assuming they 
are magnetic dipole active. For other modes the polarization dependence is not compatible 
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with magnetic or electric dipole activity of a cycloidal spin chain. We have evidence that 
there is a spin wave mode at 1 cm~^ with sufficient oscillator strength to be an electromagnon. 
Understanding the origin of the large number of optically active spin excitations and their 
interaction with magnetic field in a quantum spin system such as LiCu202 requires further 
theoretical investigation. 

We thank M. Mostovoy for useful discussions. Support by EstSF grants 5553, 6138, and 
7011 is acknowledged. Work at Rutgers was supported by the DE-FG02-07ER46382. 
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